Theory of spin accumulation and spin transfer torque in a magnetic domain wall 
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We studied the spin accumulation and spin transfer torque in a magnetic domain wall by solving 
the Boltzmann equation for spin accumulation with the diffusion approximation. We obtained 
analytical expressions of spin accumulation and spin transfer torque. Both the adiabatic and the 
non-adiabatic components of the spin transfer torque oscillate with the thickness of the domain wall. 
We showed that the oscillating component plays a dominant role in the non-adiabatic torque when 
the domain wall is thinner than the spin-flip length. We also showed that the magnitude of the 
non-adiabatic torque is inversely proportional to the thickness of the domain wall. 
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Spin-dependent electron transport in magnetic nanos- 
tructures results in many interesting phenomena such 
as the giant magnetoresistance effect [l| and current- 
induced magnetization dynamics [1, Recently, such 
spin-dependent phenomena in magnetic domain walls 
have been investigated, due to great interest in their 
potential application to spin-electronics devices such as 
spin-motive-force memory Q and racetrack memory 
[6]. In these devices, data are stored by moving the do- 
main wall using spin transfer torque. 

In 2004, Zhang and Li showed that spin transfer torque 
in a domain wall can be decomposed into two parts, the 
adiabatic and non-adiabatic torque [3]- The adiabatic 
torque lies along the spatial gradient of the local mag- 
netization while the non-adiabatic torque is perpendicu- 
lar to this direction. Assuming that the spin accumula- 
tion obeys the phenomenological diffusion equation and 
is spatially independent of the domain wall, Zhang and 
Li showed that the ratio of the magnitudes of the adi- 
abatic and non-adiabatic torques is determined by the 
precession frequency of the spin accumulation due to the 
exchange coupling and the spin-flip scattering time, and 
that the non-adiabatic torque is about two orders of mag- 
nitude smaller than the adiabatic torque. 

The thickness of a domain wall is determined by the 
competition of the exchange coupling between the local- 
ized magnetizations and the magnetic anisotropy, and is 
usually on the order of 100 nm for conventional ferro- 
magnetic metals such as Fe, Co, Ni, and their alloys. 
However, recent developments in the processing technol- 
ogy for nanostructures have allowed the production of a 
domain wall whose thickness is on the order of 1-10 nm, 
by reinforcing the shape anisotropy of a magnetic nano- 
wire [Sj] or trapping a domain wall in a current-conflned- 
path geometry [§]. These developments motivated us to 
study the transport phenomena in a thin domain wall 
[Tol [ill . [T^ . For such a thin domain wall, we cannot as- 
sume the spatial independence of the spin accumulation. 
Thus, it is important to estimate spin transfer torque by 
taking into account the spatial variation in the spin accu- 
mulation, which would be different from the estimation 
by Zhang and Li Q. Recently, Vanhaverbeke and Viret 



calculated spin transfer torque in a thin domain wall by 
numerically solving the time-dependent Larnior equation 
of the spin accumulation in a moving frame and showed 
that the non-adiabatic torque is one order of magnitude 
larger than that estimated by Zhang and Li when the 
thickness of the domain wall is comparable to the Lamor 
precession length . 

In this paper, we study spin accumulation and spin 
transfer torque in a domain wall by solving the Boltz- 
mann equation with a diffusion approximation. We ob- 
tained the analytical expressions of spin accumulation 
and spin transfer torque. Both the adiabatic and the 
non-adiabatic components of the spin transfer torque os- 
cillate with the thickness of the domain wall. We show 
that the oscillation plays a dominant role in the non- 
adiabatic torque when the domain wall thickness is less 
than the spin-flip length, which is defined by the product 
of the Fermi velocity and the spin-flip scattering time. 
For a domain wall that is much thinner than the spin- 
flip length, the non-adiabatic torque is about one order 
of magnitude smaller than the adiabatic torque, which 
is one order of magnitude larger than that estimated by 
Zhang and Li [7] and qualitatively consistent with the 
results of Vanhaverbeke and Viret [l^ . We also showed 
that the magnitude of the non-adiabatic torque is in- 
versely proportional to the thickness of the domain wall. 

We considered electron transport in a one-dimensional 
magnetic nanowire with a 180° domain wall which lies 
over —d/2 < x < d/2, where d is the thickness of the 
domain wall. We assume that the interaction between 
the conducting (s-like) electrons and the localized (d- 
like) electrons is described by an sd exchange interaction, 
Hsd — — (>^/2)o- • S, where a- is the vector of the Pauli 
matrices, J is the sd exchange coupling constant and 
S(a;) — (0, — sm9, cos 6) is the unit vector pointing along 
the direction of the localized spin angular momentum. 
The angle 9{x) is given by for x<-d/2, {-n / d){x + d/2) 
for ~d/2<x<d/2, and tt for x>d/2, respectively. 

Following Simanek and Rebel [ij) [iBl , we employ the 
rotating frame where basic unit vectors are deflned as 
Gx — a^^S X [d^/dx), By — —a~^dS/dx and = S, 
respectively. We assume that the direction of the lo- 
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calized spin varies slowly compared to the Fermi wave- 
length Ap, i.e., a = d6'/da; ^ 2t:/\y\ thus, we could 
neglect the higher-order terms of a in the following 
calculation. The spin accumulation and spin transfer 
torque in a domain wall are obtained by solving the 
Boltzmann equation for the Wigner function defined as 
f{x,Px) = [f{x,Pa;)i + g{x,Pa:)-a-]/2, whcrc f{x,px) and 
g{x,px) represent the charge and spin distribution func- 
tions, respectively. The spin accumulation s and the spin 
current density j are defined as 



d^p 
d^p 



Vxg, 
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respectively. It should be noted that the dimensions of 
s and j are density and density times velocity, respec- 
tively. The diffusion approximation, J t;^gd'^p/(27rfi)^ ~ 
(wp/3)s, is applied to the Boltzmann equation Up 
to the first order of a, the transverse components of the 
spin accumulation, Sx and Sy, and spin current, jx and 
jy, obey the following equations; 
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where luj = J/h is the Larmor precession frequency, T is 
the momentum relaxation time, Tgt is the spin-flip scat- 
tering time and D — VpT/3 is the diffusion constant 
The longitudinal spin current jz in Eq. ^ 
is given by jz = f3je/{—e), where /3 and je are the spin 
polarization factor and the electric current density, re- 
spectively. In our definition, the positive electric current 
corresponds to the electron flow along the —x direction. 

The physics behind Eqs. Q-® are as follows. Trav- 
eling through the domain wall, the conducting electrons 
vary the direction of their spin along the localized spin 
angular momentum S. Then, spin accumulation and spin 
current polarized along the y-direction (cx dS/dx) are 
induced; see Eq. The accumulated spins precess 

around S due to the sd exchange coupling with the pre- 
cession frequency ujj. Then, the x-components of the 
spin accumulation and spin current are induced [see Eq. 
([5])]. Equations ^ and @ relate the spin accumulation 
and spin current by the diffusion constant. 

Before estimating the spin accumulation and spin 
transfer torque in a domain wall, we should emphasize 
the vahdity of our calculations. Since Eqs. Q-® are 
obtained by applying the diffusion approximation to the 
Boltzmann equation, they are valid for d > Imtp, where 
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FIG. 1: The spatial variation of the transverse spin accumu- 
lations for a thick (d — 100 nm) and thin (d = 10 nm) domain 
wall; (a) Sx for d = 100 nm, (b) Sy for d = 100 nm, (c) Sx 
for d = 10 nm and (d) Sy for d = 10 nm, respectively. The 
magnitudes of Sx and Sy are divided by j^/coj, see Eq. ((Tj. 

^mfp = vfT is the mean- free-path of the conducting elec- 
trons. For a domain wall the thickness of which is much 
smaller than the mean free path, i.e., d <C ^mfp, the Boltz- 
mann equation should be solved without the diffusion ap- 
proximation. Moreover, in such a very thin domain wall, 
we cannot neglect the higher-order terms of a. 

We assume that the transverse spin accumulation and 
spin current are such that they vanish at the limit of 
\x\ — > oo and are continuous at x = ±d/2. Then, solving 
Eqs. ([SI)-®, the transverse spin accumulations in the do- 



main wall are obtained as Sx 
where 



Re[s+] and Sy = Im[s-| 
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Here C = 2/{ujjTsf) and £ is given by 
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where — Re[l/£] and fci — lm[l/£] characterize the os- 
cillation and damping of Sx and Sy due to the sd exchange 
coupling and the spin-dependent scattering, respectively. 
As shown in Eq. ([7]), the transverse spin accumulations 
can be decomposed into spatially independent (first) and 
dependent (second) parts. 

Figure [1] shows the spatial dependence of the trans- 
verse spin accumulations, Sx and Sy, for thick {d = 100 
nm) and thin (d = 10 nm) domain walls, respectively. For 
convenience, Sx and Sy are divided by jz/^j [see Eq. ([7])]. 
The parameters are taken to be J = 1.0 eV, Tsf = 10^^ 
ns, and Z^fp = 3.0 nm, respectively. The Fermi velocity is 
given by v-p — ^J2e-p/'m, where the Fermi energy is taken 
to be e-p = 5.0 eV. These are typical values for the con- 
ventional transition ferromagnetic metals . As shown 
in Figs. [1] (a) and [1] (b), for a thick domain wall, the 
spin accumulation in the domain wall is nearly spatially 
independent except at the boundaries of the domain wall 
X = ±d/2. On the other hand, as shown in Figs. [1] (c) 
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and [T] (d), for a thin domain wall, the spin accumula- 
tions vary in the domain wall, and we cannot assume the 
spatial independence of the spin accumulations. 

Let us estimate spin transfer torque in the domain wall, 
which is defined as 



d/2 



uj js X Sdx 
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The Landau-Lifshitz-Gilbert equation for the localized 
magnetization M = — S with the torque t is given by 



dM 



7M X B + aoM X 



dm 

~dt 
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TT-lTrTxm X ^— 
27rAf dx 
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where 7 is the gyromagnetic ratio, B is the effective mag- 
netic field, M is the magnitude of the magnetization and 
Q!o is the Gilbert damping constant. Ty = By ■ t and 
Tx — Sx'T correspond to the adiabatic and non-adiabatic 
torque, respectively. By using Eq. ([7]), we find that 
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The first terms of Eqs. pT|) and are identical to the 
adiabatic and non-adiabatic torque estimated by Zhang 
and Li 0, respectively. These first terms arise from the 
spatially independent part of the spin accumulation, i.e., 
the first term of Eq. ([7]). It should be noted that these 
terms are independent of the thickness of the domain wall 
d. For a thick domain wall, these first terms are dominant 
for spin transfer torque, and the ratio of the magnitude of 
the adiabatic and non-adiabatic torque, \Tx/Ty\, is given 
by C ~ 10"^ Q- On the other hand, the second and 
third terms of Eqs. pTjl and p2p arise from the spatial 
variation in the spin accumulation, i.e., the second term 
of Eq. ([7]). As shown in Figs. [T] (c) and [T] (d), for a 
thin domain wall, we cannot neglect the spatial variation 
in the transverse spin accumulation, and these second 
and third terms dominate the spin transfer torque. It 
should be noted that these terms are inversely propor- 
tional to the thickness d. Thus, for a thin domain wall, 
the strength of the spin transfer torque is considerably 
different from that estimated by Zhang and Li [t}. 
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FIG. 2: (a) The strength of the adiabatic torque Ty renor- 
malized by n/Sje/e against the thickness of the domain wall 
is shown, (b) The strength of the non-adiabatic torque Tx 
renormalized by nPje/e against the thickness of the domain 
wall is shown. For d < Zmfp = 3 nm, the diffusion approxima- 
tion cannot be applied to the Boltzmann equation, and thus, 
torque below d < Znifp denoted by the dotted line is not valid. 

Figure [2] shows the strength of the adiabatic torque Ty 
and the non-adiabatic torque t^ renormalized by 7r/3je/e 
against the thickness of the domain wall d. We denote 
the torque for d < ^mfp by the dotted line because our cal- 
culations are restricted for d > Imfp', thus, the torques for 
d < ?mfp are not valid. As shown in Fig. [21 for d > 30 nm, 
spin transfer torque is nearly independent of the thick- 
ness d. On the other hand, for d ^ 30 nm, the strength 
of the non-adiabatic torque increases as the thickness d 
decreases. For d < 10 nm, \Tx/Ty \ ~ 10"^, which is one 
order of magnitude larger than that estimated by Zhang 
and Li [7]. Moreover, for d < 10 nm, the spin transfer 
torque oscillates against the thickness d with the period 
of the oscillation given by 2-^ jki. 

Let us reveal the parameters which characterize the 
above behavior of the spin transfer torque. Assum- 
ing that k, ~ VS/iilmfp) < fci ~ VScjj/vF (Ref. [15|) 
and d ^ ^mfpi we find that Ty ~ 7r/3je/e and t^ — 
— 7rC/?je/e — 7r/3je/(ec?fci), respectively. Thus, for d^l^ip, 
the adiabatic torque is nearly independent of the thick- 
ness. On the other hand, for d<l/(Cfci)~;sf/(2V3)~40 
nm, where Isf — v-pTgf is the spin-flip length, the torque 
due to the spatial variation in the spin accumulation is 
dominant for the non-adiabatic torque. For a thin do- 
main wall, the ratio of the adiabatic and non-adiabatic 
torque is characterized by vp / {VSujjd), which is the ratio 
of the precession frequency of the electrons' spin around 
the localized spin angular momentum due to the sd ex- 
change coupling and the angular velocity of the rotation 
of the exchange field in the domain wall. For c?=10 nm, 
vp/ {VSuJjd) ~ 10^^. The oscillation period is given by 
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27r/fci ~ 2iTv-p / {s/^ito j) ~ 2.5 nm. These estimations can 
be confirmed by the plots shown in Fig. [21 

When the precession frequency of the electrons' spin 
around the exchange field, tjj, is comparable to the 
angular velocity of the rotation of the exchange field 
in space, nvp/d, the direction of the electrons' spin 
cannot vary their direction adiabatically, and the non- 
adiabaticity, which is sometimes called the mistracking 
effect, plays an important role in the spin-dependent 
transport phenomena. For example, the terminal veloc- 
ity of the domain wall motion is proportional to the ra- 
tio of the adiabatic and non-adiabatic torque The 
non-adiabaticity is characterized by a dimensionless pa- 
rameter ^ = TTVp/{2du!j) [l^. As shown above, the ratio 
of the adiabatic and non-adiabatic torques for a thin do- 
main wall, vf /{Vdujjd), is the first order of ^, while the 
magnetoresistance due to the mistracking effect or spin 
accumulation is on the second order of ^ [l^, [ij, [l^ . For 
conventional ferromagnetic metals with d > /mfp, C is less 
than unity. Thus, the non-adiabaticity plays an impor- 
tant role in the dynamics of the localized magnetization 
compared to the magnetoresistance. It should be noted 
that for a thick domain wall the non-adiabatic torque is 
characterized by C = 2/(wjTsf), not ^, as shown by Zhang 
and Li 0] . 

We compare our results with those of Vanhaverbeke 
and Viret [l3|. In Ref. [l^, a time-dependent phe- 
nomcnological Larmor equation for the magnetic moment 
in a moving frame is solved numerically, and showed that 
the non-adiabatic torque is shown to be one order of mag- 
nitude larger than that estimated by Zhang and Li [7] 
when the thickness of the domain wall is comparable to 
the Larmor precession length Al — VF/(2Trujj), which is 
on the order of a few nanometers. On the other hand, we 
consider the spin diffusion in the domain wall in a steady 
state by solving the Boltzmann equation in the rotated 
frame, and analytical expressions of the spin accumula- 



tion and spin transfer torque are obtained. We show that 
the strength of the non-adiabatic torque increases as the 
thickness of the domain wall decreases for d<lsi/{2^/3). 
Note that the condition is determined by the spin-flip 
length Isi instead of the Larmor precession length Al . We 
also find that the strength of the non-adiabatic torque 
is characterized by the first order of the non-adiabatic 
parameter ^ cx 1/d. The non-adiabatic torque does not 
change its sign, as shown in Fig. 5 in Ref. [13]. 

Since the diffusion approximation is applied to the 
Boltzmann equation, the present theory is not applica- 
ble to the ballistic region d < Imfp- The spin transfer 
torque in the ballistic region is obtained by Waintal and 
Vilet [l^. One can easily confirm that our Eq. ([7]) re- 
duces to Eq. (12) of Ref. pj] in the limit of f, Tsf oo, 
where they assume that the local spin transfer torque is 
proportional to the spin accumulation. One might ex- 
pect a simple connection formula between ballistic and 
diffusive spin transfer torque like Wexler's formula for 
conductance [13] ■ However, this is beyond the scope of 
the present paper. 

In conclusion, we studied spin transfer torque in a do- 
main wall by solving the Boltzmann equation for spin ac- 
cumulation, and found their analytical expressions. For 
a thin domain wall whose thickness is much thinner than 
the spin-fiip length, the ratio of the magnitude of the adi- 
abatic and non-adiabatic torque is about 10^^, which is 
one order of magnitude larger than that estimated in Ref. 
0] and consistent with that in Ref. [13] • We also found 
that the strength of the non-adiabatic torque is inversely 
proportional to the thickness of the domain wall. 
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